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Abstract 

The univariate spline quasi-interpolants (abbr. QIs) studied in this 
paper are approximation operators using B-spline expansions with coeffi- 
cients which are linear combinations of discrete or weighted mean values 
of the function to be approximated. When working with nonuniform par- 
titions, the main challenge is to find QIs which have both good approx- 
imation orders and uniform norms which are bounded independently of 
the given partition. Near-best QIs are obtained by minimizing an upper 
bound of the infinity norm of QIs depending on a certain number of free 
parameters, thus reducing this norm. This paper is devoted to the study 
of two families of near-best QIs of approximation order 3. 

Keywords : spline approximation, spline quasi-interpolants. 
AMS classification : 41A15, 41A35, 65D07. 

1 Introduction 

A spline quasi-interpolant (abbr. QI) of / has the general form 

where {B a ,a G A} is a family of B-splines forming a partition of unity and 
{X a ,a £ A} is a family of linear functionals which are local in the sense that 
they only use values of / in some neighbourhood of S a =supp(-B Q ). The main 
interest of QIs is that they provide good approximants of functions without 
solving any linear system of equations. In this paper, we want to study the 
following types of QIs: 

Discrete Quasi-interpolants (abbr. dQIs) : the linear functionals are linear 
combinations of values of / at some points in a neighbourhood of S Q (see e.g. 
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Integral Quasi-Intcrpolants (abbr. iQIs) : the linear functionals are linear com- 
binations of weighted mean values of / in some neighbourhood of S Q (see e.g. [2]- 

More specifically, we study Qls that we call Near-Best Quasi-Interpolants (abbr. 
NB Qls) which are defined as follows: 

1) Near-Best dQIs: assume that X a (f) = S/3gf q a a(0)f(xp) where the finite 
set of points {xp,fl £ F a } lies in some neighbourhood of S Q . Then it is clear 
that, for ||/||oo < 1 and aeA, |A a (/)| < ||a Q ||i, where a a is the vector with 
components a a (/3), from which we deduce immediately 

IIQIIoc < V \\ a {f)\B a < max |A Q (/)| < max ||a a ||i = v x {Q). 

L — ' aGA aeA 

Now, assuming that n =card(F a ) for all a, we can try to find a* £ R™ solution 
of the minimization problem 

||a* ||i = min{||a Q ||i; a a € R n , V a a a = b a } 

where the linear constraints express that Q is exact on some subspace of poly- 
nomials. Thus, we finally obtain 

IIQIloo < *f(Q) =maj Kill. 

aeA 

2) Near-Best iQIs: assume that X a (f) — J2peF a a(P) J s Mp(t)f{t)dt, where 
the B-splines Mp are normalized by J Mp — 1. Note that the B- spline Mp can 
be different from B a . Once again, for ||/||oo < 1, we have 




whence, as we obtained above for dQIs, 

IIQIloo < max||a Q ||i = v{{Q). 

aeA 

As emphasized by de Boor (see e.g. 0J, chapter XII), a QI defined on non 
uniform partitions has to be uniformly bounded independently of the partition 
in order to be interesting for applications. Therefore, the aim of this paper is 
to define some families of discrete and integral Qls satisfying this property and 
having the smallest possible norm. As in general it is difficult to minimize the 
true norm of the operator, we have chosen to solve the minimization problems 
defined above. 

The paper extends some results of £Q [TI], and is organized as follows. We 
first recall some "classical" Qls of various types and we verify that they are 
uniformly bounded. Then we define and study several families of discrete and 
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integral QIs, depending on a finite number of parameters, for which we can 
find vl(Q). We show that this problem has always a solution (in general non 
unique). Of particular interest are the results of theorems 1,4,6 and 8 where we 
show that some families of dQIs and iQIs are uniformly bounded independently 
of the partition. By imposing more constraints on the non-uniform partitions, 
we can also prove that some families of QIs are near-best (theorems 5 and 9). 
(A parallel study of spline QIs is done in for uniform partitions of the real 
line and in [3] for a uniform triangulation of the plane). In all cases, the QIs 
that we study are only exact on P2, i.e. their approximation order is 3. It 
seems surprisingly difficult to construct QIs which are both uniformly bounded 
independently of the partition and exact on Fd for d > 3. In the last section 11, 
we consider an example of QI which is exact on P3 (i.e. of approximation order 
4) and uniformly bounded: however the bound depends on the maximal mesh 
ratio of the partition. Such operators seem also useful for applications and it 
would be interesting to study near-best operators of this type, thus allowing a 
reduction of the upper bound of the norm. 



2 Notations 

We shall use classical B-splines of degree to on a bounded interval / = [a, b] 
or on I = R. For the sake of simplicity, in the case / = R, we take a strictly 
increasing sequence of knots T = {£j, 4 £ Z} satisfying |tj| — > +00 as \i\ — > +00. 
In the case / = [a, b], we take the usual sequence T of knots defined by (see e.g. 

t— m = • • • = to = a < ti < tl < ■ ■ ■ < i n -i <b = t n = • •• = t n +m- 

For J = {0, . . . ,n + m — 1}, the family of B-splines {Bj,j G J}, with support 
Sj = [tj- m , tj + i] is a basis of the space S m (I, T) of splines of degree to on the 
interval / endowed with the partition T. These B-splines form a partition of 
unity, i.e. X^eJ = 1- We set hi = U — £j_i for all indices i. 
Let N m = {1, . . . , to} and Tj = {tj_ r+ i, r s N m }: we recall that the elementary- 
symmetric functions <Ji{T) of the to variables in Tj are defined by ao(Tj) = 
and for 1 < / < m, by 

= ^ tj+l-ntj+l-ra-'-tj+l-ri- 

l<r 1 <r2<-..<ri<m 

Denoting Ci, = ,,, m! ,,. the binomial coefficients, then, for < I < to, the 

G m l\{m—l)i ' ' — — ' 

monomials e/(x) = x l can be written e; = Y^ieJ^i &h with 9^ = ai(Ti)/C l m . 

This is a direct consequence of Marsden's identity (U, chapter IX). In particular, 

the Greville points have abscissas 9i = 

Similarly, we define the extended symmetric functions <Ji(Tj) by ao{Tj) = 1 and, 
for 1 < I < to, 

= ^2 *j+X-ri*j'+X-r 2 • • • tj+i-n- 

l<ri<r2<...<ri<m 
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Then, the moments of the B-spline Mj of degree m — 2, with supp(Mj) = 
[tj-m+iftj] an d normalized by J Mj = 1, are given by (see e.g. |17j): 



3 Uniformly bounded discrete quasi-interpolants 

exact on P 2 

It is possible to derive discrete QIs from the de Boor-Fix QIs |S] by replacing 
the values of derivatives D l f(8i)/l\ by divided differences at the points 9i lying 
in £j. Doing this, we loose the property of projection on S m (I,T). However, 
by choosing conveniently the divided differences, we can obtain some families 
of dQIs which are uniformly bounded and exact on specific subspaces of poly- 
nomials. 

Let us construct for example a family of dQIs of degree m which are exact on 
P2. We start from the de Boor-Fix functionals truncated at order 2: 




W) = -Y,(- i r- l D m - l Mr)D l f(r). 
m ' — ' 



where ipj (t) = (t — t 



,+1) . . . (t — tj). We obtain successively 



D m ^(r) = (-l) m m!, D m "Vi(r) = (-l) ro m!(r - 9 3 ), 

D m - 2 ^(T) = U-l) m m\{T 2 - 29,1 + Of ). 

I J 



More specifically, taking r = 9j, we get 



D m -Vi(0 3 O - 0, D m - 2 ^{9j) = -(-l)^m\(9f - 9]) 



Thus, we can we define the QI exact on P2 



Q 2 / = ^A J (/)S ; 



whose coefficient functionals are given by 



with of ] = 9) - 9^ 



We recall the expansion (se e.g. |4][rB]1: 
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On the other hand, ^D 2 f(9j) coincides on the space P2 with the second order 
divided difference [0j_i, 9j, 9j+i]f, therefore the dQI defined by 

with coefficient functionals 

x*{f)=f(9 j )-ef ) [9 j - 1 ,9 h e j+1 ]f, 



is also exact on Pa. Moreover, one can write 



A*(/) = a j f{9 j ^) + b 3 f(6 3 ) + Cj f(9 j+1 ), with 

_i(A0 i _i+A0 7 ), bj = : 

A9j). So, according to the introduction 



-0f ) /A0 J _ 1 (A0 J _ 1 +A0 3 ), b = 1+9^ /A9 j . 1 A9 j , Cj = -9f /A9 j (A9 j ^ 1 - 



IIQ2II00 < max(|a,-| + \bj\ + |cd) < 1 + 2 max 6»f ) /A9 j ^ 1 A9 j . 
je.j jeJ J 

The following theorem extends a result given for quadratic splines in lll | |22 | P3 ) . 

Theorem 1 For any degree m, the dQIs Q\ are uniformly bounded. More 
specifically, for all partitions of I : 

||05lloo<[i(m + 4)]. 

proof: We only give the proof for m = 2k + 1, the case m = 2k being similar. 
For the sake of simplicity, we take j — m, i.e. we shall determine an upper 
bound of the ratio 

N m /D m = 9<$/A9 rn _ l A9 rn 

with 

N = 0( 2 ) = - V (t - t) 2 

m z {m — 1 ' 

v ' \<r<s<m 

Setting H = J^_ 2 then we get a lower bound for the denominator 

1 1 H 2 

Dm = -^(t m - t )(t m+1 -h) = ^(hi+H)(H+ h m+1 ) > —. 
m z m z m A 

The numerator N m is composed of k pairs of sums (S pi S' p ) 

S p = Yl (tr-tsf, S' p = (tr-t*) 2 , 
s — r—p s — r—m—p 

for 1 < p < k. Both sums contain at most p times the terms hf and 2hihj 
(i 7^ j), hence we can write S p + S' p < 2pH 2 , which implies 



N m < 



2H 2 {k + l)H 2 



-(l + 2 + ... + fc) = ., 

m z (m — 1) Zm z 
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so, we get 

N m /D m < i(fc + l), 
and finally, for m = 2k + 1 odd 

HQ2II00 < k + 2 = him + 3) = [him + 4)]. 

For m — 2k, we obtain respectively D m > ^ and N m < 4 ^_ 1 - ) , whence 
N m /D m < 2 ^_ 1 , and finally for m = 2fc even 

||Q^||oo<fc + 2=i(m + 4) = [i(m + 4)]. □ 

4 Existence and characterization of near-best dis- 
crete quasi- inter polants 

4.1 Existence of near-best dQIs 

We consider the following family of dQIs defined on I = M endowed with an 
arbitrary non- uniform strictly increasing sequence of knots T = {tj, i G Z}, 

Qf = Q P , g f = J2 Xi ^ Bi - 

Their coefficient functional depend on 2p + 1 parameters, with 2p > to, 

p 

M/) = £ *(«)/(*<+«). 

and they are exact on the space P g , where q < m. The latter condition is 
equivalent to Qe r = e r for all monomials of degrees < r < q. It implies that 
for all indices i, the parameters ai(s) satisfy the system of q+1 linear equations: 

s— — p 

The Vandermonde matrix Vi £ H>(9+ 1 ) x ( 2 p+ 1 ) of this system, with coefficients 
Vi(r, s) = 9 r i+s , is of maximal rank q + 1, therefore there are 2p — q free parame- 
ters. Denoting bi € R 9+1 the vector with components bi(r) — 9^ r \ < r < q, 
and by dj G M 2p+1 the vector with components aj(s), we consider the sequence 
of minimization problems, for i e Z: 

min||aj||i, V i a i = b i . (Mj) 

We have already seen in the introduction that ^i(Q) = maxig^niin ||cii||i is an 
upper bound of ||Q||oo which is easier to evaluate than the true norm of the dQI. 
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Theorem 2 The above minimization problems (Mi) have always solutions, which, 
in general, are non unique. 

proof: The objective function being convex and the domains being affine sub- 
spaces, these classical optimization problems have always solutions, in general 
non unique. □. 

4.2 Characterization of optimal solutions 

For b € K m and A S W mxn , let us consider the Zi-minimization problem 

min ||r(a)||i, r(a) = b — Aa. (M) 

We recall the characterization of optimal solutions for ^-problems given in |3(J|. 
chapter 6. Define the sets 

Z{a) = {l<i<m; n{a) = 0}, 

V(a) = {v£ W n ; IHU <l,Vi = sgn{n(a)) for i { Z(a)}. 

Theorem 3 The vector a* G K™ is a solution of (M) if and only if there exists 
a vector v* £ V(a*) satisfying A T v* = 0. 

5 A family of spline discrete quasi-interpolants 

exact on P2 

In this section, we restrict our study to the subfamily of spline dQIs which are 
exact on P2, i.e. we consider the dQIs Q p = Q p ,2- 

We shall need some set of indices. Let K = {— p, . . . ,p}, and K* = {— p, 0,p}. 
Then, we can write K := K \ K* — K\ U K2, where K\ = {— p + 1, . . . , —1}, 
and K2 = {1, . . -,p— 1}. 

The three equations expressing the exactness of Q p on P2 can be written 
a i (~P) + a i (0) + a% (p) = 1 - ^ a-i (r) 

r£K 

0i- P a l (-p) + B i a l (0) + 8 l+p a,i (p) =6j— 8 i+r ai (r) 

reK 

el pai (-p) + efa (o) + e* +p a t (p) = - ^ (r) 

Let (a* (— p) , a* (0) , a* (p)) be the unique solution of the system with the right- 
hand side obtained by taking a, (r) = for all reK. Using Cramer's rule, we 
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obtain 



a* (-p) = -ef } / (6 i+p - 9^ p ) (9i - 6i_ p ) , 
a* (0) = 1 + ef } / {9 l+p - 6i) (6i - Oi-p) , 
a* (P) = -e? ] l (0 i+P - 0i- P ) (0 i+P - 0t) ■ 
Then we can express the general solution of the above system in the form 

a, (-p) = a* (-p) - a ^ a i ( r ) + ai ' sCli ( s ) ' 
di (0) = a* (0) - Y /V ffl i ( r ) _ X! A,« a i ( s ) ' 
a 4 (p) = a* (p) + Y 7»,rO» (f) ~ X! 7i, s ai (s) , 



with 



and 



A. 



(Oj — Qj+j ) (9j+p — dj+j ) 

{@i — @i-p) {Oi+p — 6i-p) ' 

(Oj+j - Pi) (Oj+p - Oi+j) 

ifii — Oi-p) (0i+ p — Oi-p) 

(Oj+j — Qj-p) (Oj — Oj+j) 

[Oi+p ~ ®i-p) {Oi+p ~ Qi) 

(Oi+j — 0i- p ) (0 i+ j — 6i) 
{Qi+p — Oi-p) {0i+ p — 0i) ' 

(Oi+j — 0i- p ) (0i +p — Qi+j) 



(0i — Oi-p) (Oi+p — 9i) 
We denote by Q* the spline dQI whose coefficient functional are 

A* (/) = a* (-p) f (0^) + a* (0) / (0,) + a* (p) f (0 i+p ) . 

In that case, an upper bound of the norm of this QI is max ie z v* , where 

vt = \a*(-p)\ + \a*(0)\ + \a*(p)\. 

Theorem 4 For all p > m, the infinity norms of the spline dQIs Q* are uni- 
formly bounded by ^±1 . This bound is independent of p and of the sequence of 
knots T. 

proof: We have to find a good upper bound of u*. We recall that 0i = 



SreN m U+i-r and 0i 



(2) 



m 2 (m— 1) 



Si, with 



Si — (ti+i- r — ti+i- 8 ) 

Kr<s<m 
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Define the following sums: 

m 

S% = (U+l-r — U+l-m-r) = hj+2-r-k, 

r£N m rGN m fe=l 

m 

(ti+l+m-r — ii+l-r) = fe-i+l-r+fc- 
r£N m reN m fc=l 

As p > m, we obtain 



1 1 p 1 

#i — = — / , (U+l-r — U+l-p-r) = — / , / , ^-i+2-r-fc > — 5*2, 

m ^ m ^— ' to 

1 1 P 1 

— @i = — / (U+l+p-r — U+l-r) = — / / , hi+l+k-r > — S3. 



m to ' '-^ to 

rENm r6N m fc=l 

The proof being essentially the same for all i e Z, we can restrict our study to 
the case i = to. In that case, we get 

m— 1 



5 2 = m/ii + fc (hm+i-k + hk+i- m ) > S' 2 = (to - 1) ft 2 H h 2/i TO _i + /i„ 

m— 1 

5*3 = mh m+1 + ^2 k ( h 2m+i-k + h k +i) > S' 3 = (m — 1) h m -\ h 2/i 3 + ft 2 . 



fc=i 



Setting Hk = h,2 + ■ ■ ■ + hk+i, for 1 < fc < m — 1, and H = H m -\ for short as 
in the proof of theorem 1, we have 



i=0 

whence 

S2S3 > S' 2 S' 3 = mH ^ - I H i } ■ 

i=0 V i=0 / 

Now, we come back to S\ and we shall prove that Si < S' 2 S' 3 < S2S3. Si can 
be written under the form 

m—lm—r m — 1 m — 2m — 1 

5 * = E E + • • • + = E + E E - H rf - 

from which we deduce 

m—l m — 2 m— 1 m—1 /m— 1 \ ^ 

5 1 = (m-i)^^-2E^ E ^ = *»E fl ?- E ff < • 

i=l j — 1 i=j+l i—1 \ i= 1 / 



Then we use the fact that, for all 1 < i < m - 1, Hi < H, whence Hf < H H t 



and YlT=o H? < Hj^iLo Hi- we obtain 



i=0 



. i=0 



Finally, for all i G Z, we have 



"•* = 1 + ;^-n& sl 

whence ||Q*|| < maxi/* < m ~*~ \ □ 



m + 1 



m — 1 m — 1 ' 



In the next section we prove that the quasi-interpolants Q* are near-best in the 
sense of section 4 under some additional conditions on the partitions. 



6 The family Q* of discrete quasi-interpolants is 
near-best 

Let us write the minimization problem (Pd) of section 4 in Watson's form. 
Taking into account the expression of the solution ai of the system equivalent 
to the exactness on P 2 of Q p 2, we can write 



Mli = Ik* - AidiW , 



where 



a-i = (a { (-p + 1) , . . . , a { (-1) , a* (1) , . . . , a* (p - 1)) G 



a* - (a* (-p) , 0, . . . , 0, a* (0) , 0, . . . , a* (p)) 1 G 
and Ai G R(2p+i)x(2 P -2) ig giyen by 



a i,-p+l 






• a i;P _i 


-1 














-1 






















-1 














-1 


-7i,-p+l ' ' 


• -7i,-i 


7i,i 


Ji.p— 1 



Theorem 5 Assume that the sequence of knots T satisfies, for all ieZ, the 
following properties 



-1 + @i < #i-p + #p < + 
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then, for all i 6 Z, a ( * is an optimal solution of the local minimization problem 
{Mi). Thus, for all p > m, the spline dQIs Q* (theorem 4-) are near-best. 

proof: According to theorem 3, we must find a vector v* £ M. 2p+1 satisfying 

Halloo < 1, Ajv*=0, v* (r) = sgn(a* (r)) for r = -p,0,p. 

Let us choose v* (— p) = — 1, v* (0) = 1, v* (p) = — 1, and 

-ai t r + Pi,r + 7»,r, H j £ Kl, 
Cti,r + 0i,r - 7t,r, if 3 £ K 2 . 



Then it is easy to verify that the equations Af v* — are satisfied. Moreover, the 

—(2) 

above expressions of a* (r) for r = —p,0,p, with i > imply that v* (r) = 
sgn(a* (r)). It only remains to prove that \v* (j)| < 1 for j £ K\ U i^2- As 
A,r = 1 - °V + 7»,r f° r r e .Ki and /3 M = 1 + a i;S - 7», s for s e if 2 , it is 
equivalent to prove 

< a i>r - 7 i>r < 1, < 7 i;S - a i;S < 1, for (r, s) £ K\ x if 2 . 

We only detail the proof for r £ K\, that for s £ Ki being quite similar. Using 
the expressions of on^ r and 7^ r given in section 5, we get 

(0% - Oi+r) [{Oi+p + 9i-p) - {Oi+r + Si)} 



&i,r 7 



l.r 



(6i — 6i- p ) {Oi+p — 0i) 
and we shall have — 7i, r > if and only if 

Oi+r + 6i 5: $i+p + $i-p 

for all r £ K\. However, since we have 9i+ r + 6i < + there only remains 
the unique condition 

+ &i < $i-p + #j+p- 

The other inequality a^ r — 7^ r < 1 can be written 

(9i — 9i+ r ) [{9i+ p + 6i- p ) — {9i+ r + 9i)\ < (9i — Qi- P ) (0i+ p — 6i) . 

Setting 5\ = 9i+ r — 9i- p , 82 = 9i — 9i+ r , and ^3 = 9i+ p — 9i, the latter inequality 
can be written 82 {83 — 81) < 83 {82 + Si), or equivalcntly 81 (82 + £3) > which 
is obviously satisfied. For s G K2, the inequalities < 7i, s — cti, s < 1 are satisfied 
if and only if 

Oi-p + 9i+ p < 9i + Oi+i, 
whence the conditions on the sequence of knots. □ 

Remark 1 Theorem 5 imposes some additional conditions on the sequence of 
knots. For quadratic splines, we have studied arithmetic and geometric se- 
quences: in both cases, the higher is p, the stronger are the conditions and, 
for p — > +oo, T is closer and closer to a uniform sequence. 

Remark 2 Even if the partition T does not satisfy the hypotheses of theorem 5, 
the operator Q* is still a good QI because its infinity norm is small and uniformly 
bounded. 
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7 Uniformly bounded integral quasi-interpolants 
of Goodman- Sharma type 

General integral spline quasi-interpolants (iQIs) already appear in Jl] |2J • 
Here we have chosen to study a family of QIs which we call Goodman-Sharma 
(GS-) type iQIs, as they first appear in |l(Jj . They seem simpler and more in- 
teresting than those studied in |2J and • I n the case of splines of degree to 
on a partition in n subintervals of a bounded interval / = [a, b], the simplest 
GS-type iQI can be written as follows: 

n+m— 2 

Glf = f(t )B + M/)£i + /(*»)£»+m-l, 

i=l 

where the integral coefficient functionals are defined by 
Ai(/)= / Mi(t)f(t)dt = (Mi,f), 



Mi being the B-spline of degree to — 2 with support Sj = [ti- m +i, U], normalized 
by A, (e ) = /Lt- = J R Mi = 1. It is easy to verify that G\ is exact on Pi and 
that HGilloo = 1- In this section, we shall study the family of GS-type iQIs 
defined by 



ra+m— 2 



Ghf = I (to) B + faA-i (/) + (/) + CiXt+i (/)] Bi+f(t n ) B n+m _i, 

»=i 

which we want to be exact on P2. The three constraints G2&k = efc, k = 0, 1, 2, 
lead to the following system of equations, for each l<i<n + m — 2: 

(2) (2) (21 (2) 

a* + bi + a = 1, 0i_iaj + 6»j6j + 0j+iCj = 0,, /i-J 1 a l + fi\ 'b t + ^ + i Ci = °l ■ 

We recall the values of the first moments of Ma: 



1 ST 1 sr 



H-i — / L i+l-r — / H—m+r 

to * — ' to - — ; 

l<r<m l<r<m 

(2) 2 v -v 2 v -v 

Mi = 7 ; 77 / ^ ^j+l— rti+l— 8 = 7 j 77 / H-m+r'i-m+8' 
TO ( TO + 1 ) TO ( TO + 1 ) 

v ' l<r<s<m y ' l<r<s<m 

Theorem 6 For any degree m, the iQIs G2 are uniformly bounded indepently 
of the partition of I . For to = 2k or 2k + 1, there holds 

IIGalL <2fc + 3. 

proof: Taking the differences of the second and third equations above (Gi&k = 
efc, k = 1, 2) with the first one (G^eo = eo) times and /ij resp., we get 



(2) /)(2) 
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Setting di = (9i+i — 9i)cn and Cj = (6i — ai, we obtain 

A/i^Aft - Apf } A0i_i] ai = ^ - 0< 2) . (E) 

Using the expressions of the various coefficients in terms of symmetric functions 
of the knots, we obtain the following form for the coefficient of Qj in equation (E) 

2 

[(*»+! - U+l-m) Act 2 (Ti-l) - (U - ti- m ) Act 2 . 



m 2 (m + 1) 



Setting l i =t i - i i+ i_ TO , we can write i f - = ft i+ i_ m + - ij 

fti+i-m + k + h i+ i and - i J+1 _ m = + hj+i- Then we have 



i+l 



A(7 2 — — ti + i_ m ) ^ t r , 



r— i — m+1 

and the coefficient of in equation (E) is given by 

AMi-lA#i — A^ A#j_i = — (U+l — ti + l- m ) (ti — ti- m ) {U+l — ti-m) 

= — (h + h i+ i) (/ij_|_i_ TO + Zj) (/ij + i_ m + + /ij+i) • 
Now, writting a 2 (T-) = er 2 + X^=i+i-m*r' we get successively 

A4 2) - ^ = -^-^^ (T t ) - -r^a* (TO 
' " 1 m (to — 1) 

((to -1)ct 2 (Ti)-(m + l)CT 2 
((m-1) ^ t?-2tf2(T0j 



TO 


(to - 


-1) 




2 




TO 


(to 2 


-1) 




2 




TO 


(to 2 


-1) 




2 




TO 


(to 2 


-1) 



m<r<s<i 



Setting Wj = Ei+i-m<r< s <i ~ ^) 2 ' we obtain 

TO 



a,- 



to 1 (hi+i-m + k) (hi + i- m + k + hi+i) (k + hi+i) ' 
from which we deduce 

1 LOi 



TO - 1 (hi + i- m + h) (hi + i- m + h + h i+ i) ' 
1 OJi 



to - 1 (h i+ i- m + k + h l+1 ) (li + h i+ i) ' 

and 

bi = 1 - ai - Cj, |a,| + |&i| + |cj| = 1 + 2 (|aj| + |cj|) . 



13 



On the other hand, we have 

1 u>i (/ij + i_ TO + 2li + hi + i) 



\ai\ + Ci = 



m-1 (ft i+ i_ m + k) (hi + i- m +k + h i+ i) (k + h i+ i) 



As LOi < k(k + 1) if for m = 2k + 1 (sec proof of theorem 1 with respect the 
upper bound for N m ), we obtain 

uji < k (k + 1) (h i+ i- m + h) (h + h i+x ) . 

Moreover, it is obvious that 

hi+l-m + 2li + hi + \ < 2 (hi+l-m + k + hi+l) . 

Therefore, we finally obtain for all i 

2k(k + l) , _ 1, 
kl + |c«|< J t _ 1 ' =fe + l = -(m + l), 

whence the uniform upper bound for the norm of G2 

l|G 2 || oc <m + 2, 

which is both independent of the partition and of the (odd) degree of the spline. 
For m = 2k even, a similar computation leads to the uniform bound 

l|G 2 || oc <m + 3. 

□ 



8 Existence and characterization of near-best in- 
tegral quasi-interpolants 

Now, we consider the family of iQIs 

whose coefficient functional depend on 2p + 1 parameters: 

p r 
A, (/) =J2 ai / Mi + S (*) / (*) d< ' 

and which are exact on P q . As in section 7, Mj is the B-spline of degree m — 2, 
with support Sj = [^ + i_ m ,^-] normalized by J R Mj = 1. The constraints 
(e r ) = e r are equivalent to the following systems of linear equations, for all 
i eZ: 

2 /x£.U(«) = *i r) . 0<r< 9 , 
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whose matrix coefficients Wi £ r(9+ 1 ) x ( 2 p+ 1 ) ; defined by W t (r, s) = fi^2 a , is of 
maximal rank q + 1 (see e.g. 26 |), thus there remains 2p — q free parameters. In 
view of the introduction, as max.; e z \\di\\i is an upper bound of the true infinity 
norm of the iQI, we want to solve the minimization problems, for all i: 

minllajll-L, Wia, = 6j. (M,) 

As in section 4.1, the objective function being convex and the domains being 
affinc subspaces, we can conclude 

Theorem 7 The above minimization problems (Mi) have always solutions, in 
general non unique. 

As in section 4.2, we shall use the characterizacion of optimal solutions for 
Zi-problems (M) given in [30], chapter 6. 

9 A family of integral spline quasi-interpolants 

exact on P2 

In this section, we restrict our study to the subfamily q — 2 of the above spline 
iQIs which are exact on P2. Moreover, we assume that p > to in order to insure 
that the three sets of knots Ti_ p , X^, and Ti +p are pairwise disjoint. Now, the 
matrix coefficients of the linear system equivalent to the exactness of G Pi 2 on 
P2 is of maximal rank 3, and we have 2p — 2 free parameters. Let us denote by 
G* the spline iQI whose coefficient functionals are 

A* (/) = a* (-p) (Mi_ p , /} + a* (0) (M t J) + a* (p) (M i+P , /) , 

where 

a* (-p) + a* (0) + a* (p) 
6i- P a*i {-P) + ^a* (0) + 6 l+p a* (p) 
$} p a* {-p) + nf ] a* (0) + n [ S v a* (p) 

that is, their coefficients are the unique solution of the system obtained by taking 
ai (r) = for all r e K (we use the same notations as in section 5). 
The following result is the analog of theorem 4 for iQIs. 

Theorem 8 For any degree m > 2, and for all p > to, the iQIs G* are uni- 
formly bounded independently of the partition of I . More specifically, there holds 




for to even, 
for to odd. 



1- 



3 (2) 
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proof: We shall prove that \a* (-p)\ + \a* (0) + \a* (p)\ < 1 + \C (m) for all 
igZ, which is sufficient to insure the result. For the sake of simplicity, we can 
assume that i — m. Solving the corresponding linear system, we get 

a m (-P) = T . a m (p) = , a (0) = l-a m (-p)-a m (p) , 

where 

f — // 2 ) - and S — A A // 2 ) - A # A // 2) 

sm — P'm u m i ana u m — L-\pU m —pL-^p^L m L-^pU m L-\pfj, rn _p : 

with 

\0i = 9 l+p - 6 U A p ^| 2) = - M; (2) , l = m,m-p. 
As proved in theorem 6, we have 

K ' K ' l<r<s<rn 

For the expression of <5 m , we need some additional notations. For 1 < i < m, 
we define 



r=2 



Similarly, for 1 < j < m, we define 

^m+j — t-m+p ~~ tj, l m +j = tj+P ~ = lm+j + ^ ' 



m — 1 



(2") 

Taking into account the definitions of 9i and ^\ for i = m — p, m, m + p, we 
obtain after some algebraic calculations the following expressions: 



_^ m _^ m 

A p r „_p = — h> A p TO = — E^ ^m+ji 
i=l j=l 

z=l \r— z s— 1 / 

7 = 1 \ r=? s=l 



Now we can write 



J = l 4=1 
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where D m is equal to 



m / m 



E^J lE'U E^+Ew - E'U EmE^+E'- 

l j-l yr— j 6—1 / / \ 7 — 1 / \i— 1 \ r— i s— 1 

mm I rn j m i 

= x X ^m+j I X] tr + X] ~ X * r -p ~ X * s 

i—l j—1 \ r—j 5=1 r= i s— 1 

Therefore, as we shall see later (p. 18) than D rn > 0, we obtain 



i=l 



where 



'2m \ /2m \ 

E'H w ™= (Em E (*r-*.) 2 - 

vfe=l / \fe=l / l<r<s<m 



Let us compute an upper bound for N m . Let H = Y]"L 2 as m ^ ne P r0 °f °f 
theorem 1. In the first sum, we have + l', m+i = k + h+ m + H , for 1 < i < m, 
so we can write 

2m 2m 

E4 = E^+ mff ' 

fc=i fe=i 

For uj m , we have already seen (also in the proof of theorem 1) that 

iv m < c(m) H 2 , 

where c (m) = k 2 for m = 2k and c (m) = k (k + 1) for in = 2k + 1. So, we 
finally obtain the following upper bound for N m : 

N m < c{m)H 2 (^h + mH^j . 

Now, we will compute a lower bound for D m . Let 

m j m « 

r— j s— 1 r—i s—1 

be the coefficient of ■ in the double sum defining D m . Therefore we have 

m— 1 m 
D m > ^m^m+l-^m,m+l + l i l Jn+1 Li^ n+ i + l m l m+ jL m ,m+j ■ 

i=l j=2 

We first observe that 

I'ml'm+lLm^+l = (lm + H) (l m +l + H) (l m + i TO+ i + H ) > 2H 2 (l m + l m+ i)+H 3 . 
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Then, we obtain successively, for 1 < % < m — 1, 

m m 

— ^ ^ tr ~t~ ^ ^ t r —p ^ (^m ^m— 1— p) ~t~ (^p+1 ^m— p) ^ 2J^. 

r— i+l r— i 

Similarly, for 2 < j < m, 

j i-i 

Lm.m+j = ^ ^ ^p+s ~~ t"m—p ~~ ^ ^ ^ (^p+1 ~~ ^m—p) (^p+2 ~~ ^l) ^ 2i?. 
s— 1 r— 1 

On the other hand, we also have successively 

i i 

KC+i = (h + E M(Wi + h) > H(h + $>), 

r=2 r=2 

m m 

cc +j = (i m + ff )(Wj+ E h s)>H{i m+] + e h ')- 

s=j + l s=j+l 

From these inequalities, we deduce 

(2m m— 1 i m— 1 m— 1 

E^ + EE^ + E E h » 
/c=l i=2 r=2 j=2 s=j+l 

Now, it is easy to see that 

m— 1 i m— 1 m— 1 

EE^ + E E >». = (m-2)ir, 

i=2 r=2 j=2 s=j+l 

therefore we obtain the lower bound 

/2m \ 

D m >2H 2 \ Yj k + (m- -)H \ . 

Thus, setting C — J2k™i 'fc' we nave the two inequalities 

I r i I ™ i u i j\r ^ „^™\ u2 



D m > 2i^ ( C + [ m - - ) H ) , N m <c(m)H z (C + mH) , 



from which we deduce 



N m 1 £ + m# 

< -c(mj- 



For m > 2 even, it is easy to verify that m(m — 1) < (m + 2)(m — 3/2), whence 



£ + mff m + 2 



£+(ra-§)i2" _ m-l 
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For to > 3 odd, one can verify that m(m — 1) < (to + l)(m — 3/2), whence 

£ + mH m + 1 



£+(m-|)i?-m-l' 
Finally, as ||G*|| is bounded above by 

2 /V 

l<4 (-p)I + l<C (0)1 + K (p)| = i + 2 (K (-p)I + Ka (p)\) < i + z— T7T ' 

(the same upper bound is valid for \a* (—p)\ + |o* (0)| + \a* (p)\ , for all i e Z), 
we obtain respectively 

^ ( m + 2)c(m) , 

G* < 1 + ' \ ' for to even 

(m - 1) 

^ 1 , ( TO + l)c(m) 
G* < 1 + „ for to odd 

ii p ii — / 1 \2 

(TO — 1) 

As c(m) = |to 2 for m even and c(m) = j(m 2 — 1) for to odd, we obtain 

\\G* p \\<l + \c{m) 

to 2 (to + 2) (to+1) 2 . 

with G(m) = — — — for to even and C(m) = — for to odd, which 

(m — l) 2 to - 1 

proves the theorem. □ 



10 The family G* of integral quasi-interpolants 
is near-best 

We follow the notations and techniques used in section 5 for discrete QIs. As 
the linear system satisfied by the coefficients of Aj (/) 

p 



s= — p 

is of maximal rank, it can be written as follows 



a i (~P) + a i (0) + a i (P) = 1 - J! ai M ' 

Oi-pai {~p) + Sitti (0) + l+p ai (p) = Oj - ^ 9i +r a t (r) , 
A4 V (-P) + M< 2) * (0) + ( P ) = 0< 2 } - J2 to 
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and its general solution is 

a, (-p) = a* (-p) - a i,r a i ( r ) + X] a ^ sCbi ( s ) ' 

a i (0) = a* (0) - X jli^ai (r) - X /3 iiS Oj (s) , 
a i (p) = a* (p) + X Ji, r ai (r) - ^ 7»,sa» ( s ) i 

where a* (— p), a* (0), and a* (p) are the coefficients of G* studied in section 9 
above, and the various coefficients are quotiens of determinants. Specifically, 

a r = W (r, 0,p) /W, p r = W {-p, r,p) /W, lr = W (-p, r, 0) /W, 
a s =W(0,s,p)/W, f3 s = W(-p,s,p)/W, >y a =W(-p,0,8)/W, 

where W (k, I, m), fc < I < m, is the determinant with columns ^1, 6* i+ fe, Mi+fc) ; 

^1, 0j + ;,/4+;) , and ^1, #i +TO , Mi+m) > an d W = W (—p, 0,p). As in section 5, 
we can write the minimization problem P; described in section 8 (with q = 2) 
in Watson's form, so we have to minimize ||ai || x = \\a* — Aicii^, where we have 
used the same notations as in section 5. 

According to theorem 3, we must find a vector v* G K 2p+1 satisfying 
IMI <1, Ajv*=0, v*(r) =sgn(a*(r)) forr = -p,0,p. 



Let us choose v* (— p) = — 1, v* (0) = 1, v* (p) = — 1, and 
v* (j) = 



-a i>r + (3i, r +1i,r, if j & Ki, 
Oti, r + Pi,r - Ji,r, if j S K 2 . 

Equations Afv* = are satisfied. Moreover, as proved in theorem 6, /i| 2 ' ) — 
6f^ > holds, and the explicit expressions for a* (—p), a* (0), and a* (p), similar 
to those obtained for i = m in the proof of the theorem 8, imply that v* (r) = 
sgn (a* (r)). It only remains to prove that, for (r, s) G K\ x 

l w * ( r )\ = \-®i,r + Pi.r + li,r\ < 1, K («) I = |ai, s + - 7i,s I < 1- 

As /3 iir = 1 - on^ r + ji >r for r e Ki and /3 i;S = 1 + a i;S - j iiS for s G K 2 , it is 
equivalent to prove 

< ai,r - 7»,r < Ij < 7 M - a i;S < 1, for (r, s)eifiX if 2 . 

We only detail the proof for « = m and for r G ifi, that for s G K 2 being quite 
similar. 

We prove that 7 m>r < a TO , r , r G ifi by stating that 7 m ,- r < a TOj _ r , for r G 
{1, . . . ,p — 1}, the latter being equivalent to W (— p, — r, 0) < W (— r, 0,p). By 
expanding the various determinants involved, we obtain the inequality 

(2) _ (2) (2) _ (2) 

l^m—r ^ Mm+p Mm-p 
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For 1 < i < m, let w> = £j =i+1 _ r h h w ' = YJjLi w 'f T l = (E^-r *i + £}=i *;) , 

With these notations, we can write succesively 

1 m 1 

m ' m 

i— 1 

1 m j 

m m 

m 

(2) _ (2) _ Z \ ^ / / 

Mm M ro -r- m(m + 1) Z^^ T i' 

o m 
(2) _ (2) _ f 

i— 1 

and the inequality 7 m ,-r < «m,-r is equivalent to 

^ m ^ m 

i— 1 i— 1 

It can be interpreted as follows: the barycenter of the m points t[ with weights 
w\ is less than or equal to the barycenter of the to points r, with weights Wi. 
The inequality a m - r — J m ,-r < 1, »" € {1, . . . ,p — 1} is equivalent to W (— r, 0,p) — 
W (—p, —r, 0) < W and can be written 

(2) _ (2) (2) _ (2) 

r^m—r r^m—p ^ r"m+p r*"m—r 

7n — r — m —p OjnJf-p — m — r 

Using the same techniques as above, we obtain successively 



^ m 

E{ti— r ti—p ) 
... 

2=1 

_^ m 

— £ fe+P - *i-r) 



J m—r v m —p 

to ^— ' m 

1=1 



1-/ 

— 10 , 



1 _ 

—w, 



where 



with 



7 m+p u m — r 

to 1 — ' TO 

i=l 



w' = w-, and w = ^ 

i=l i=l 



i+r i+p 

= /ij, and w 4 = fy, 

j=i+\-p j=i+l- r 
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for 1 < i < m. Defining r\ = ^_ (g£l p tj + ^%\_ r t 3 ) , and t> 
^TT (l)£=l-r h + Ej=i +P tj) for 1 < * < m, we get 



J 2 ) _ u ( 2 ) - t VWr' 

,(2) _ u ( 2 ) = ? V 
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Then the inequality a m ,_ r — J m> - r < 1 can be interpreted in a geometric form 
as follows: the barycenter of the m points with weights «J- is less than or 
equal to the barycenter of the m points Tj with weights ?Uj. 
In a similar way, one can prove that inequalities < ~f m ,s — c*m,s < 1, for all 
1 < s < p — 1 are equivalent to the following inequalities: 

(2) _ (2) (2) _ (2) (2) _ (2) (2) _ (2) 

^m+p t^m-p ^ Hm+s I 1 ™ f-m+s Mm-p ^ t^m+p Pm+8 

@m+p @m—p @m+s @m @m+s @m—p @m+p @m+s 

and can also be interpreted in terms of barycenters of knots. Here we need the 
weighted points (r/', w'{ ) and (t", w"), where r/' = ^ {j2™ =l tj + Ejti+s *i) ' 

r'i = ^i (e?=Zs t> + y:;=\ + p tj) , < - r^ +1+s and < = v 

Theorem 9 Assume that the sequence of knots T satisfies, for all i G Z, the 
following properties: 

1. for all 1 < r < p — 1, the barycenter of the m points (r t ',io^) (resp. 
(r-,uJ-)J is less than or equal to the barycenter of the m points (r^w^) 
(resp. {TuWi)). 

2. for all I < s < p—1, the barycenter of the m points (r i7 Wi) (resp. (fi,Wi)) 
is less than or equal to the barycenter of the m points (t" ,w") (resp. 

Then, for all i G Z, a* is an optimal solution of the local minimization prob- 
lem (Mi). Thus, for all p> m, the spline iQIs G* of theorem 9 are near-best. 

Remark 3 Even if the partition T does not satisfy the hypothesis of theorem 9, 
the operator G* is still a good iQI because its infinity norm is uniformly bounded 
(theorem 8). 



11 Quasi-Interpolants exact on P n with n > 3 

While we succeeded above in characterizing some families of near-best QIs exact 
on P 2 whose norms are uniformly bounded independently of the partition, it is 
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surprisingly difficult to find QIs exact on P3 having the same property. Actually, 
we did not find any example of such a QI. 

For example, let us consider the following cubic spline dQI (which also appears 
in |2(J|). defined on / = R endowed with a non- uniform partition, by 

Qaf(x) = ^Ai(/)B,-(af), 

where Bi is the cubic B-spline with support [ti-2,U+2\ centered at ti and 

h{f) = a>ifi-i + hfi + Cifi+i, 
where /, = f(ti), and the coefficients are given by 

1 ftf , _ 1 (hi-i + Kf 1 hti 



3 hi-i(hi-i + hi) 3 hi-ihi 3hi(hi-i + hi) 

It is easy to verify that Q3 is exact on P3, i.e. that ai,6j,Cj satisfy the system 
of linear equations: 

a,- + 6* + Cj = 1, ij-ia, + t l b l + tj+ic, = 0;, 
This operator can also be written in the quasi-Lagrange form 

iGZ 

where the fundamental function Bi, having support [ti— 3, ii+a], is defined by 

Bi = Ci-iBi-i + biBi + di+iBi+i. 

The Chebyshev norm \A\00 of the associated Lebesgue function A = J^iez 1-^*1 
satisfies: 

|A|oo = HQ3II00- 

For x G I = [£2,^3], we have A = J2i=o 1-^*1 an< ^ we shall now construct a 
partition for which \A\<x, is arbitrary large. We choose ft 3 = h as parameter and 
hi = 1 for alH 7^ 3. Then the fundamental functions on the interval I are given 
by: 

1-41-14 h 2 

Bo = - — B^ . — —Bi — —Bo. Bo — — B^ -\ — Bo — — : — B%, 

6 ' 3 6 6 3 3(1 + ft) 

1 (1 + ft) 2 1 1 
- 1 l(ft + l) 2 = 1ft 2 

Bi -~3h(l + h) B3+ 3^^ B4 > B5 -'30Th~) Bi - 
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Let us compute the value of A at the midpoint s = \{t 2 + t 3 ) of /. 

Setting ai = Si(*i) = B 4 {t 3 ), a 2 = B 2 {h) = B 3 (t 3 ), 6 2 = B 2 (t 3 ) = B 3 {t 2 ) and 

using the algebraic properties of B-splines, we get 

h 2ft 3 + ft 2 + 9 . ft 

a l = o/i , 7 T 1 "2 = , 7T779 1 , 02 



3(1 + ft)' 3(ft+l)(/i 2 -ft + 3)' " (ft + l)(ft 2 -ft + 3)' 

(wc have ai + a 2 + ^2 = 1 because the B-splines sum to one). Now, using the 
values of the central BB-coefhcients (3 2 and 72 of B 2 on the interval I: 

72 = ft 2 -ft + 3 ' A = l-T4, 
we can compute the values of B-splines at this point 

B 1 (s) = lai = B 4 («), B 2 (s) = B 3 (s) = \{a 2 + 3/3 2 + 3 72 + 5 2 ) = 1(4 - ai). 

o 00 

After some algebraic calculations, we obtain the asymptotic behaviour of the 
Lebesgue function at the midpoint of I : 

A(s)=0(ft), h->+oo, 

therefore the norm of the associated QI is unbounded. 

However, if we now assume that there exists r > such that the partition 
satisfies 

1 , fti + l , • _ 

- < — : — < r, l G Z, 
r ftj 

then we obtain the following upper bounds: 

\a t \,\ Cl \<±j^- y N<l(l + 2 , 
from which we deduce 

\\Q 3 \\oo<N(r):=U(l + r) 2 + . 2 "" 



3 V (1 + 0/ 

For example, for r = 1, 2, 3, 4, 5, we get the following values of the upper bound 
of the norm : 

N(l) « 1.66, JV(2) w 3.89, AT(3) « 6.83, JV(4) w 10.47, JV(5) w 14.78, 

which are still of reasonable size. Therefore it seems that such QIs are interesting 
in practice though they are not uniformly bounded with respect to all partitions. 
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